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Abstract 

We compute the time delay of strong lensing in the framework of the Einstein-Straus solu- 
tion. We compare the theory to the observational bound on the time delay of the lens SDSS 
J1004+4112. 
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1 Introduction 

Before the work of Rindler and Ishak [1], the general believe was that the deflexion angle 
of light passing near an isolated static and spherically symmetric mass is independent of the 
cosmological constant. This believe was based on the argument that the cosmological constant 
disappears from the geodesic equation for massless particles. In September 2007, Rindler and 
Ishak [IJ corrected this believe. They pointed out that it is not sufficient to consider the geodesic 
equation but also the metric itself must be considered. In the work of Rindler and Ishak the 
source emitting the light and the Earth were supposed to be at rest with respect to the lens. 
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Also all the masses, including those of the Earth and of the source, were neglected, except the 
mass of the lens. Since then, there is a rich controversy about whether or not a cosmological 
constant modifies the bending of light near an isolated spherical mass. Sereno [2l[3], Schiicker 
[11[5], Miraghaei and Nouri-Zonoz [6j, Kantowski, Chen and Dai [7], confirm Rindler and Ishak's 
result, while Khriplovich and Pomeransky [8J, Park [9], Gibbons, Warnick and Werner |10] . 
Simpson, Peacock and Heavens [llj contradict Rindler and Ishak's findings. 

Recently, Schiicker |12) redid the calculations of the bending of light by a spherically sym- 
metric mass distribution, which is taken to be a cluster of galaxies by relaxing all the previous 
hypotheses except that of sphericity: the observer and the source are allowed to move with re- 
spect to the cluster, the masses of the other clusters are included in the form of a homogeneous 
isotropic dust and the observer as well as the source are taken comoving with respect to the 
dust. The appropriate framework of this computation is the Einstein-Straus solution \1'6\ [Hj 
that matches the Kottler solution at the inside of the Schiicking radius with the Friedmann 
solution at the outside. Einstein and Straus' first motivation was to explain why the cosmic 
expansion does not affect smaller length scales such as planetary and atomic systems. Schiicker's 
calculations confirm previous calculations by Ishak, Rindler, Dossett, Moldenhauer and Allison 
|15j : Taking into account realistic cosmic velocities attenuates the effect of the cosmological 
constant on the bending of light without however cancelling it. Although there has never been 
a claim that the time delay was independent of a cosmological constant, it is also interesting to 
do the computation of the time delay in the framework of the Einstein-Straus solution. 

We will use the same units as in reference [T^]: astrometers (am), astroseconds (as), astro- 
grams (ag) 

am = 1.30 • 10^6 m = 4221 Mpc, 

as = 4.34 • 10^7 s = 13.8 Gyr, (1) 
ag = 6.99- 10^1 kg = 3.52 • 10^^ M©, 

where Mq denotes one solar mass. In these units 

c = lam as~^, SvrG = 1 am^as^^as^""^, Hq = 1 as^^ , fi = 3.86 • 10^"'^^"'^ am^ as^^'^ag. (2) 

For spatially fiat universes, to which we will restrict ourselves in the following, we may set the 
scale factor today oq = 1. 

2 The Einstein- Straus solution with a cosmological constant 

We will consider hereafter the Einstein-Straus solution \\.'6\ [H] generalized to include the cosmo- 
logical constant [16] but restrict ourselves to spatially flat universes. We will need the Jacobian 
of the transformation passing between the Friedmann and Schwarzschild coordinates to calculate 
the geodesies of photons. Let us quote the results obtained by Schiicker [12]. Let (T, r, 6, if) and 
{t, X, 0, if) stand respectively for Kottler and Friedmann coordinates. The Kottler metric 

ds^ = B{r)dT^ - B{r)-^dr^ - r^dQ^, (3) 

with 

2GM A p 

B{r) = 1 -r^ 4 

r 3 

prevails inside a vacuole of radius rschuiT), r < rschu- The Friedmann spatially flat metric is 
given by: 

ds'^ = dt^ - a{tf (dx^ + x^dif) , (5) 
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with the scale factor a{t) determined by the first order Priedmann equation 



da 
'dt 



A/a + Aa^/3, 



where 



^ = O'OPdusto/'i = 1 



A 

3' 



(6) 
(7) 



prevails outside the vacuole x ^ XSchu- It is worthwhile to notice that due to ([6]), the scale 
factor is strictly monotonic. The two solutions are glued together at the constant Schiicking 
radius xschu 

rschu{T) ■■= a{t)xschu- (8) 

By taking into account the fact that the central mass must be equal to the dust density times 
the volume of the ball with Schiicking radius rschu 



A -- 

B{rschu) = 
It is also useful to introduce Cschu 



2M 



2GM 



B 



A 



Schii 



XSchu 



A 



2 2 
" XSchil- 



Cschu '■= V 1 - Bs, 



chii' 



(9) 
(10) 

(11) 



Schiicker [I2j computed the Jacobian of the coordinate transformation at the Schiicking radius 
and also the inverse of the Jacobian (corresponding to the Jacobian of the inverse coordinate 
transformation {t,x) ~^ {T,r)) with the results 



and 



dt_ 

dT 

dT 

dT 
'dt 

dr 
'di 



Schii 



Schii 



Schii 



Schii 



C Schii 
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B Schii ' 
C Schii ) 
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Cschii 
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Schii 


Bschii 


dx 
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dr 


Schii 


aB Schii ' 


dT 




Cschii 

a 

Bschii 


dx 


Schii 


dr 
dx 


Schii 


a. 



(13) 



In the following we also need to pass between Kottler time T and Friedmann time t at the 
Schiicking radius. To this end we will use the result also obtained by Schiicker [T2] 



dt 
dT 



Schii 



Bs, 



chii- 



(14) 



3 Integrating the geodesies of light 

We have the following situation: a first photon is emitted by the source, a quasar, at a time t'g 
and follows an upper straight line trajectory until its arrival at a time fg^f^^^, on the Schiicking 
sphere in the half space containing the source. It is then bent inside the Schiicking radius, until 
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it emerges from the Schiicking sphere in the half space containing the Earth at time t^c/iiiE' then 
follows a straight line until its arrival on Earth at time = 0. A second photon is emitted by the 
quasar at a time ts follows a lower straight line trajectory, arrives on the Schiicking sphere in the 
half space containing the source at time tschuSi it is then bent inside the Schiicking sphere and 
emerges from the Schiicking sphere in the half space containing the Earth at time tschUE, follows 
again a straight line until its receipt on Earth at the same time as the first photon tQ = I'q = 0. 
We will here be interested in the computation of the time delay ts — t'g. 




Figure 1: Two light rays emitted by a source S, bent inside the Schiicking sphere and finally 
received at Earth E. The travel times of the two photons differ, giving rise to a time delay. 



Let us first integrate the first order Friedmann solution ([6]) for the scale factor a{t) in the 
spatially flat case with cosmological constant A and dust density Pdusto = (3 — A) (we will take 
the experimentally favored value A = 0.77 • 3 am~^ it 20%) with final condition a(0) = 1. One 
may show that ^ is equivalent to the Friedmann second order equation for a{t) 



2 (fa ]_f^Y_. 



with final conditions 



da 
IE 



t=o 



a(0) = 1. 



(15) 



(16) 



The integration of ([6]) or ()15p may be done analytically with the result 



ait) 



1 - A/3 y/^ 



sinh 




We will also need to solve the equation 



dx 
dt 



t + arcsinh 



1 

a 



A/3 
1 - A/3, 



1/3 



(17) 



(18) 



for various final conditions, with x having the meaning of a geodesic distance, not to be confused 
with a luminosity distance. Since these are the final conditions, at the arrival on Earth, which 
are known, we will proceed backwards in time in three steps: we will determinate t'g^y^^^ and 
tschuE, then t^chiis and tschiis, and finally t'g and ts- 
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step 1: Determination of t'g^f^^^ and tgchuE- 

Here we will be interested in the photon trajectory between the Earth and the Schiicking 
sphere with the Friedmann metric. The non vanishing Christoffel symbols of the Friedmann 
metric in the plane 6 = 7^/2 are given by 

= -X, r,^x = (19) 

rf^ = at/a, r^^ = l/x, 

and the geodesic equation reads 

i+aatid + aatx^ip^ = 0, 



a 

+ 2— + -xv? = 0, 
a X 

with final conditions at p = 

* = 0, X = XB, (y2 = TT 



(20) 



i = 1, X = cos a', 



sin a' 



XE 



(21) 



for the upper trajectory photon, where we use the fact that the physical angle a' coincides 

with the coordinate angle arctan(|x0/xl)- 

For the lower trajectory photon, the final conditions differ 



t = 0, X = XE, if = -TT 



sma 

t = 1, X = cos a, If = 



(22) 



XE 



The solution of the geodesic equation is 



i = ^ = sin(<^-a'), 'P = (23) 

a X « X 

where x'p, the would be peri-lens, is given by 

Xp = XfiSino;', (24) 

for the upper trajectory photon and 

i = ^ = -sin(<^ + a), ^ = (25) 

a X ax 
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where Xpj the would be peri-lens, is given by 

Xp = XEsma, (26) 

for the lower trajectory photon. 

The polar angles ^'schuE VSchuE at which the lower and the upper trajectory photons 
emerge from the Schiicking sphere are given respectively by 

f'schuE = - arcsin ( — ^ ) + a' , (27) 
\ XSchil J 

fSchiiE = -vr + arcsin ( ] - a, (28) 



XSchil 



where XSchu, the Schiicking radius, is given by 



/My/3 

^'^'"^ = Uvr(l-A/3)j ' 

where we make use of d?]) together with the fact that we work in system of units such that 
A = 1 — A/3. It is easy to see using elementary Euclidean geometry that the geodesic dis- 
tance x'schiiE E between the upper trajectory photon when it emerges from the Schiicking 
sphere and the Earth is given by 



XSchuE,E = \/Xe + Xlchu + 2 XEXSchu COS f'schilE^ (30) 

where xe is the Earth-lens geodesic distance. 

In the same way the geodesic distance XSchuE,E, between the lower trajectory photon when 
it emerges from the Schiicking sphere and the Earth is given by an analogous expression 



XSchuE,E = V X| + Xschil + 2 XEXSchu COS ifSchuE, (31) 

The Earth-Lens geodesic distance xe may be deduced from the knowledge of the redshift 
ze according to the scheme 

ZE ^ aE = > Te ■■= t{aE) x(tE) =■ XE, (32) 

i + Ze 

where t{a) denotes the inverse of the scale factor a{t), where x is the solution of the 
first order differential equation p8p . subject to the initial condition x(0) = 0, meaning 
that the photon reaches the Earth today at to = 0. Moreover, since the scale factor a is 
strictly positif, X is a strictly decreasing function and thus injective. Therefore, x{t) can 
be inverted to give t in term of x- 

The knowledge of ze allows to determine qe, then one deduces the corresponding time tE- 
Injecting in x(^)) one finally deduces xe- This is the same procedure that will be used to 
deduce xs 

zs ^ as = — > ts ■■= t{as) x(ts) =■ Xs- (33) 

I + 

On the other hand, the geodesic distance between the upper trajectory photon and the 
Earth in the time interval from the crossing of the Schiicking sphere in the half space 



6 



containing the Earth t'g^y^^^ until its arrival at Earth = is governed by equation (fT8|) 
with the final condition 

X(0) = (34) 

meaning that at tg = the photon reaches the Earth. From (|18p and (j34p . one deduces 
that 

''^SchilE.E f^SchuE 1 

dx = - -77\d^ (35) 

JO a[t) 

i.e., 

, ft'schilE 1 , /"O 1 , , , 

XschuE,E = - / = / ^dt. (36) 

•^0 a{t) Jt',^,^^a{t) 



Then, comparing with equation (j30p . one gets 

V X| + xLhii + 2 Xi?XSch« cos (/j's^/^^^ = —-dt. (37) 

Equation ([37|) may then be used to deduce t'g^^^^. 

The geodesic distance between the lower trajectory photon and the Earth in the time 
interval from the crossing of the Schiicking sphere in the half space containing Earth, 
tschuEi until its arrival at Earth (to = i'o = 0), is also governed by (fT8|) with the final 
condition (fM]l . In an analogous way, one gets 

/•o 1 /-o 1 

XSchuE,E = / -nrdt = / -7-:dt. (38) 

■JtschilE a(rj -JtsohuE «(*) 

Then, comparison with equation (j3ip gives 



/-O 1 

x| + x|c/i« + 2 XEXSchu COS <y95ch^i£; = / -TTTf^i, (39) 



which may serve to deduce tschuE- 

However, we find it more reliable to proceed in a different way: we compute tschuE by 
difference with t'g^^^^. Combining (f37|) and ([39]) . one gets 

, /-o 1 , /-o 1 /-^'scMB 1 , 

XSchuE,E - XSchuE,E = / ^TTOt - / ^TTfi^t = / -JT^dt. (40) 



tscMB a{t) Jt's.^^j, a{t) JtschuE a{t) 

Since a{t) does vary significantly only on cosmological time scales and since |t5(,h«_E; ~ ^SchHEl 
is very much smaller than cosmological time scales then 

tSchuE '^ytschuE) 
On the other hand, combining (|30p and (|3ip 



XSchuE.E - x'schuE,E = V x| + Xschil + 2Xi?XSchii COS (/?Sch«£; (42) 

Xe + Xschil + '^XEXSchu COS ^'^c;,^^;. (43) 
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But 



cos tpschiiE = - cos ( - arcsm — ) + a 

\ \Schu 



cos 1 h a 

XSchu 

g - X. 

XSchu 



and similarly 



, , 1 / XE - XSchu\ 2 /^rX 

COS LP schuE ^ + ^\ a, (45) 

and then one gets up to second order in the physical angles a and a' 

XSchuE,E - x'schiXE,E - iXE - XSchu) „ ("^ - a'^) ■ (46) 

^XSchu ^ ' 

Combining (|^T|) and ([Ml), one gets an approximate expression for t'g^^^^^ — tschuE 

t'schiiE - tschiiE ^ a{t'schuE)iXE - XSchii) ^ (a^ - a'^) , (47) 

which may be used to deduce tschuE if ^'schuE been determined from ([57]) . Since 
a > a', the lower trajectory photon emerges from the Schiicking sphere before the upper 
trajectory photon. 

The upper trajectory photon emerges from the Schiicking sphere with 4-velocity 

II _ 1 ■/ _ cos {^'schuE - «0 ■/ /.ox 

^SchilE — ~ 1 XSchuE — a ' ^SchuE — ~P2 2 ' 

"'SchilE "'SchuE "'SchuE^Schu 

where 

^'schuE '■= ^i't'schiiE)^ (49) 

and the lower trajectory photon emerges from the Schiicking sphere with 4-velocity 

; 1 . COS {if SchuE + a) . -Xp 

tSchuE — , XSchilE — 2 ' ^SchuE — — 2 ' V-^U) 

aSchuE (^SchuE '^SchuEXschil 

where 

aSchilE '■= a{tschuE)- (51) 

Let j'p, 'jF be the smaller physical angles between the un-oriented direction of the upper 
trajectory photon and the direction towards the lens, and between the un-oriented direction 
of the lower trajectory photon and the direction towards the lens. We have 



7p = arctan 

and 

= arctan 



XSchu ■ I 

^SchilE 

^SchilE 



XSchil' 



XSchuE 



- {^'schuE - "0 = arcsin I I (52) 

V XSchii I 



) = TT + {(pschiiE + a) = arcsin ( ) . (53) 

/ V XSchu / 
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step 2: Determination of i^c/iuS' ^^'^ tschuS- 

We have first to translate the 4- velocities of the upper and lower trajectories photons into 
the coordinates (T, r, (p) . They will serve as final conditions for the geodesic equation inside 
the Schiicking sphere where prevails Kottler metric. Using the inverse Jacobian ()13p . one 
gets 

^> _ C'schuE - cos {^'schiiE - 



' SchiiE 



and 



with 
and 

where Cschu{t) is given by 
Cschu{t) = 



^SchilE 

CschuE - COS {ipSchilE + a) 
aSchilE 

^SchiiE = Cschuit'schuE) 

C SchiiE = CschilitSchuE), 



1 — Bschu{t) 



(54) 

(55) 

(56) 
(57) 

(58) 
(59) 



Let 7^ and denote respectively the smaller coordinate angles between the un-oriented 
direction of the upper trajectory photon and the direction towards the lens and between 
the un-oriented direction of the lower trajectory photon and the direction towards the lens. 
We have 

sin 7^ 



and 



7^ =: arctan 



7x ='■ arctan 



„/ ^SchuE 
' SchiiE ■/ 

' SchilE 

'f SchiiE 



r SchiiE - 



rSchiiE 



arctan 



arctan 



C'schiiE + ^OS-f'p 
C SchiiE + COS 7F 



(60) 



(61) 



Moreover, we have at our disposal an initial condition, [12], which we may use to set 
'^SchiiE ~ ^SchiiE- have thus specified the final conditions of the geodesic equation 
inside the Schiicking sphere. 

Making use of the Christoffel symbols of the Kottler metric in the equatorial plane 9 = Tr/2 



BB' 



TT 



2B' 



-rB, 



-TV 
^ rip 



the geodesic equations then read 



r + -B{r)B'{r)t'^ 



1 B'{r) 

2 B(r) 



•• B'(r) ■ 
B[r) 

— rB{r) ip^ = 0, 



ip + -ftp = 0, 
r 



(62) 

(63) 
(64) 
(65) 
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from which we deduce three first integrals 



T=l/S(r), (66) 
^r'^ = J, (67) 

+ -T-^.=-E. (68) 



B{r) B[r^ 

Equation ()67p comes from invariance of the metric under rotations and J has the meaning 
of an angular momentum per unit mass, ()68p comes from invariance of the metric under 
time translations and E has the meaning of energy per unit mass. For the photon E = 0. 

Eliminating the affine parameter between (j67p and (|68p . one gets in the case of a photon 



At the peri-lens for the upper trajectory photon ^ , = 0, from which one deduces the 
expression of J in terms of r'p 



Similarly at the peri-lens for the lower trajectory photon ^ 
another expression for J in terms of 



and one deduces 



J= . (71) 



Birpj 

Replacing J by the appropriate expression ([7D|) or ([7T1) . one gets 

^ = ± i f 1 - ^ - (72) 

valid for the upper trajectory photon and 



a(p 1 / 1 •s \ 



(73) 



valid for the lower trajectory photon, where s denotes the Schwarzschild radius s = 2GM. 

It is worthwhile to notice that the cosmological constant has disappeared from (|72p and 
also from (1731). 



Since s/r'^ <C 1 and s/rp <C 1, we will hereafter only retain terms up to linear order in 
s/r'p or in s/vp. In this approximation 

r'p ~ a{t's^j^i^E)xschu sin 7^ - GM (74) 

and 

rp =i a{tschuE)xschu sin 7/^ - GM. (75) 
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Eliminating now the affine parameter between (|66p and ()68p and taking into account ()70p , 
one gets 

■ ^ (76) 



B{r) 



for the upper trajectory photon and 

dT _ 
dr 



BirJl 



(77) 



r2 B{rp) 



for the lower trajectory photon. 

Let us denote T^c/iiiS '^'schuE Kottler times at which the upper trajectory photon 
penetrates inside and emerges from the vacuole respectively and ^5^/1^5 ^'schUE 
corresponding Friedmann times. From now on, we will denote the expression of dip/dr 
and dT/dr valid for the upper trajectory photon by dip' /dr and dT' /dr respectively and 
will continu to denote the expressions valid for the lower trajectory by d(p/dr and dT/dr. 
We have 



SchilE ~ ^ SchuS 



t 

P 



dT' 



dr 



dr + 



dT' 



dr 



dr, 



(78) 



with r{t) = a{t)xschu and where dT' /dr is given by ([76]) . To obtain ([78|) . we have used the 
fact that T' decreases when r increases from r{t'g^f^^g) to r'p as well as when r increases 
from r'p to rit's^^^^). 

Using the relation (|14p relating the Kottler time T and the Friedmann time t, one gets 



'^'schuS — '^'schuE 



dt 



Bschilit) 



(79) 



Moreover, we have at our disposal a free initial condition [12] which we may use to set 
^'schiiE ~ '^SchilE instance. Combining ([78]) and ([79]) . one gets 





dT' 


dr + 




dT' 


dr — 




dr 




dr 


Ir'p 


u 


Jt' 



SchilE 



dt 



Bschiiit) 



0, 



(80) 



from which one can deduce t'^chuS- interested in T'g^j^^^ one can use ([7^ to obtain 

it. 

Let us denote TschuS and TschuE the Kottler times at which the lower trajectory photon 
penetrates into and leaves the Schiicking sphere respectively and by tschuS and tschuE the 
corresponding Friedmann times. 

To compute tschuSi one can proceed in a similar manner as for t'^^^^^g. We have the 
analogous of ([TH]), (I79D and 



TschilE — TschuS 
TschuS = TschuE 



r{tschuE) 
rp 

tSchiiE 



dT 


dr + 




dT 


dr 




dr 





dr, 



dt 



tSchilS 



Bschu{t) 



(81) 
(82) 
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and 



r{tschuE) 
rp 



dT 



dr 



dr + 



r{tschus) 
rp 



dT 



dr 



dr 



dt 



0, 



(83) 



from which one can deduce tschiis- 

But as in the cas of tschuE, we prefer to proceed in a different manner: determining tschuS 
by differences. Combining (|78p and (|8ip. one obtains 



i'^'schuE ~ '^'schus) 





dT' 


dr — 




dT 








Jr' 

V 


dr 


Jrp 


dr 



+ 

On the other hand, using ([H]), we have 



{TschuE — Tschus) - 
dr 



^ SchuE' 



ritschiis) 



dT' 



dr 
dT 



dr 



dr 
dr. 



(84) 



and 



'^SchiiE ~ TschuE 



'^SchiiS ~ TschilS 



SchilE 



dt 



'^'schuE ^SchilE 
Bschu{t'schuE) 



tSchilE Bschiiit) 

dt 

tschuS Bschu{t) Bschu{t'schus) ' 



i'c 



SchiiS 



tSchilS 



(85) 



(86) 



where we have used the fact that Bschu does not vary appreciably on the time intervals 
[tschuEit'schuE^ i^SchuS ^t'schus\ since these are smaller than cosmological scales. 
Substituting ^ and ^ into 



^'schuE tSchilE 
Bschuii'schuE) 



(|Mj) one gets 




^'schuS tSchuS 


Bschu{t'schus) 


r'^^'^'schuE) 


dT' 




dr 


rr{tschuE) 


dT 




dr 


f'^i'^'schuE^ 


dT' 


Jr'p 


dr 


r^^y'-'schuS^ 
+ / 


dT' 

dr 


rr(tschuE) 


dT 


^ SchuE' 


dr 



dr + 
dr — 

dr — 
dr 



I 

V 

ritschus) 
P 

^'(^'schilE') 
P 

P 

r{tschus) 
r(t' ,..„) 



dT' 



dr 
dT 



dr 



dT 



dr 
dT 



dr 
dT 



dr 



dr 
dr 

dr 
dr 
dr. 



(87) 



But since we deal with smaller length and time scales than cosmological ones: 

[a{tschuE) - a(*!s'c/in£:)] XSchil, 



rr{tschuE) 


dT 


dr ~ 


dT 






^ SchuE'' 


dr 


dr 


rr{tschus) 


dT 


dr ~ 


dT 




Jr(t' , .-e) 


dr 


dr 



^ bchuE' 



r(t' . ..„) 

SchuS-' 



[a{tschus) - a{t'schus)] XSchu- 
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Using equation ([6]), one deduces 
a{tschuE) - a{t'schuE) - 
a{tschus) - a{t'schus) ~ 

with 



A A 

+ IT^^'^i^'schuE) i^SchuE - t'schuE) i 



^i'^'schuE) 



A 



a{t[ 



SchiiS' 



+ -^"-^i^'schiis) {iSchuS - t'schus) 



A = - 



^PdustfiCLQ — -pdustfi- 

On the other hand, using the results of reference [12j 





dT' 


dr — 




dT 


M 


"l 


Jr'p 


dr 


o 1 p 


dr 


^ Svr 


2 



.5 a-i 

2 I r{? 



M 



87rr2^]V73 



It suffices to make the fohowing replacements 



87rr2 




arctanh J i7r{t'schuE) " 21n M 



(90) 
(91) 



(92) 



X — )• — , eT — i 

in equation (23) of reference [T^. 
In the same manner 



''"('^'schilE) 



M 




87rri 



21n(^ 



(93) 



(94) 



Collecting together the previous results, one gets an analytical approximation for tgchuS ' 



SchilS 



P 



3 I r: 



^-1 



M 



arctanh [ J ^r{t'g^f^^j^) j + arctanh 



2 yr'p^ 



-41n ^ 



' '^^('^'schus) 



dT 
dr 



^^(t'scHnE) iXE - XScnn) - «'^) 



dT 




dr 


'■(t'schiis)) 



(95) 



13 



where we have also used the expression (|T7|1 for {t'g^f^;^^ — tschuE)- Then the knowledge of 
t'schus allows one to deduce tschuS- 

Let us now determine in turn the polar angles ^'schuS ™^ ^SchuS at which the upper 
and lower trajectories photons penetrate inside the Schiicking sphere. Since the angle (p' 
increases all the way from r'^^^^^^ to and from to ^'^cftiiE' S^ts 



V'SchuS — ^'schuE 





dip' 


dr — 




d^' 


Ir' 

P 


dr 




dr 



dr. 



(96) 



To the linear order in the ratio s/r^, Schwarzschild radius s = 2GM divided by peri-lens 
r'p, one gets [15] 



V^SchilS 



VschuE - ^ + arcsin 



SchuE , 



+ arcsm 



' SchuS , 



' SchuE 



Ir'c 



SchuE 



r. 
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' SchiiE 



' SchuE 



2 r' 



l_s_ 

2r' 



r'schilS 

r.l2 



r' 

SchilS 



SchilS 



+ r' 



(97) 



Let us now compute the polar angle fschuS at which the lower trajectory photon penetrates 
inside the Schiicking sphere. Since ip decreases when r varies from rsckuS to Vp and also 
when r varies from to rschuE^ then 



rr(tschus) 


dip 


dr + 




dip 


dr. 




dr 


Ir-p 


dr 





^SchilS = ^SchilE + 



To the linear order in the ratio s/rp, one gets for pschuS [E 



V^Schiis - V'SchiiE + TT - arcsin 
M 



rp \ 



TSchilE , 



arcsm 



rp \ 



r SchilS , 



+ 



1 



Svr 



rschiiE \ 



' SchuE 



1 + 



M 



Svr rschuS \ 



SchilS 



M 1 



/ rschilE 



M 1 r SchilS 



rschiiE + rp Svr rp V rschiis + rp 



(98) 



(99) 



step 3: Determination of t'g and tg 
One can now compute ip'g 



"Ps = ^SchilS - iFS + arcsin 



/ XSchil . I 
sm7^5 



(100) 



where ip'schns given by ([97|) and where the smaller physical angle between the un- 
oriented direction of the photon and the direction towards the lens as the photon penetrates 
inside the Schiicking sphere from the external side: 



i'fs 



arctan 



^ '^SchilS 
XSchil ■ f 

XschilS 



arctan \-XSchii 



^'schilS \ 
XschilS I 



(101) 
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where ^p'schuS ^nd x'^chuS are given respectively by 

r' 

p 



^'schus = V= = -— 7^ 7==, (102) 



and 



x'schilS 



Cschujt'schus) \ 



(103) 



with 



and 



/ A A 

Cschuit'schus) = J^T^T—^Xschu + ■^'^Ht'schus)Xschu (104) 

/ / \ 2 2 / / \ 2 / \ 

Bschuitschiis) = 1 777 TXSc/iu ~ "T" V' SchuSlXSchu- (105) 
"'V'SchilS) 

To obtain ()103p we have used the Jacobian of the coordinate transformation (T, r) — )■ (t, x)) 
(fTB|) . together with the expressions of T^c/iuS ^^"^ '^'schuS 



1 



T'schuS — 1T777 V' ''^SchiiS 



'2 

SchilS) 



(106) 



Xl_5 is the geodesic distance between the source and the lens, which can be accurately 
approximated by 

XL,S ~XS - XL- (107) 

In the same manner, once (pschuS determined, one can compute the polar angle ips corre- 
sponding to the source by a relation analogous to (jlOOp 

fs = VSchiis + IFS - arcsin gin jps ] , (108) 

V XL,S J 



where ^pschuS is given by (j99p and where ^ps is the smaller physical angle between the 
un-oriented direction of the lower trajectory photon and the direction towards the lens as 
the photon penetrates inside the Schiicking sphere from the external side: 



IFS = arctan 

Using 

XSchilS 



SchilS 

XSchii ■ 

XSchilS 



) = arctan ( XSchu^^^^^) ■ (109) 

/ V XSchilS J 



^ _ Bschiisjt Sch ilS) 

0.{tschils)B SchilS itSchils)\ '''SchilS B{rp) 

C SchilS jtSchils) 
aitschils)BschilsitSchils) ' 



(110) 
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obtained by using the Jacobian of the coordinate transformation (T, r) — t- (i, x) and 



'fSchilS = ; (111) 



'{tSchus)Xschu^jBirp) 

then 



7i?S = arctan 



'^pBschusi'tSchus) 



, rschuS\/B{rp) i CschilsitSchus) + \h 

\ V V ^Schii 



(112) 



For a given M, one obtains in turn 1'^^,^^,^, t^^hiiS' <^Sc/i«5' tschiiE, tschiis, VSchiis and 

In general 93'^ 7^ fs- To achieve 93'^ = which corresponds to the fact that the upper 
and the lower trajectories photons are emitted by the same source, we have to adjust M, 
i.e., we have to vary M until the equality ip'g = ips is satisfied. We end up with values of 

^'schuE^ t'schuS^ f'schilS^ tSchiiE, tSchiiS, fSchilS and tp'g = ips- 

We are now in a position to determine ts — t'g. 

Using once again some elementary Euclidean geometry, similar to that used to obtain 
x'sch'iiE E (1301) > one obtains for the geodesic distance x'schuS s between the source S and 
the photon of the upper trajectory as it crosses the Schiicking sphere in the half space 
containing the source: 



x's,SchuS = ^Xl,s + Xschil - '^XL.SXSchu COS {^'schus " fs)- (113) 

Proceeding in the same manner, we get for the geodesic distance XSchuS,s between the 
source S and the photon of the lower trajectory as it crosses the Schiicking sphere in the 
half space containing the source 



XSchuS,S = yXL,S + Xschii - '^XL,SXSchu COS {(fSchilS " Vs)- (114) 

Making use of the approximations 

cosx~l — 2;^/2 and \/l + x ~ 1 + x/2, (115) 
valid for |x| <C 1, one gets 

Xschus,s - iXL,s - Xschu) i^-zi TTT' '-'^^") 

2 \Xl,S ~ Xschu ) 

I X i}PSchuS - ^sf /, -, „N 

XSchuS,S - \XL,S-XSchu) ITT' ^ > 

'^\XL,S-Xschu) 

Then 

XschuS,S - XSchuS,S —T ■ (ll»j 

[Xl,S - Xschu) 
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On the other hand, equation (fT8]l with the final condition x{t'schus) ~ *-* gives 

0-x'schuS,S = -l -7TT, (119) 

Jt'g a[t) 

In a similar manner, from (jlSp with the final condition xitschus) = 0, one gets 

ftSchilS 

0-XSchuS,s = - (120) 
Jts a(t) 

From (|119p and ()120p . one deduces an expression for xschuS,s — x'schuS s 

^ rtschiis dt r^'schiis dt 

XSchiXS,3 — XSchuS,S = / "TTT ~ / " TTT 

Jts 0,{t) Jt'g a{t) 



*'s dt j-tschiis dt r^'schiis dt 
ts «(0 Jt' a{t) Jt' a{t) 



s ^ ^ "s 

s 



dt ftschiis dt 



But since a(t) varies significantly only on time intervals of cosmological nature 

*S jt_ ^ t'g-ts rtsMS ^ tschiiS - t'schuS ^22) 

, a{t) <^KcH^s a{t) a{t's,f,^s) ' 

Then 

xschus,s - Xschus,s ~ —77rr 7777 n — • (i^-s) 

a[Zg) "'V'Schiis) 
Equating the right hand sides of (jllSp and (jl23p . one arrives at 

- t's _^ tschuS - t'schuS ^ {y^SchilS - Vsf - {'P'schuS ~ 'Psf ^-^24) 

Then one deduces an expression for ts — t'g 
ts - t's ~ ait's) 



t'schilS - tSchuS ^ i'PSchuS - ¥^5) - if'schilS ~ 'PsY 



<i'schus) 2 (xl]s - Xslhu 



(125) 



Hereafter are displayed the results for ips{= f's) the deflexion angle, M the mass of the 
cluster of galaxies (the lens) and At = ts — t'g, the time delay, in the case of the lensed quasar 
SDSS J1004+4112 where [Ml [IH [20]. '±0' stands for the central value, '+' and '-' stand 
respectively for the upper and the lower experimental limits. 

a' = 5" ±10%, Q = 10" ±10% . . 

ZL = 0.68, zs = 1.734. ^ 

The cluster mass M comes from a fitting: the angles ips and ip'g are calculated as a function of 
M which is varied until the equality <p'g = ips is satisfied. 



17 



For the cosmological A, we take the experimentally favored A = 0.77-3 am ^ it 20%. We have 
also considered the case without cosmological constant A = in table HI 



Table 1: Upper limit value of A: A = 0.77 • 3 am'^ + 20% 



aE±10% 


±0 


±0 


±0 


+ 


+ 


+ 








a'e ± 10% 


±0 


+ 




±0 


+ 




±0 


+ 






9.03 


8.13 


9.94 


10.84 


9.94 


11.74 


7.23 


6.32 


8.13 


M [10^3^0] 


1.80 


1.98 


1.62 


1.98 


2.18 


1.78 


1.62 


1.78 


1.46 


[years] 


9.76 


9.18 


10.25 


12.35 


11.81 


12.77 


7.38 


6.74 


7.91 



Table 2: Central value of A: A = 0.77 • 3 am"^ 



ae ± 10% 


±0 


±0 


±0 


+ 


+ 


+ 








a'e ± 10% 


±0 


+ 




±0 


+ 




±0 


+ 






9.97 


8.98 


10.97 


11.97 


10.97 


12.97 


7.98 


6.98 


8.98 


M [IO^^Mq] 


1.82 


2.00 


1.64 


2.00 


2.21 


1.80 


1.64 


1.80 


1.48 


At [years] 


9.72 


9.14 


10.19 


12.28 


11.76 


12.68 


7.35 


6.73 


7.87 



Table 3: Lower limit value of A: A = 0.77 • 3 am~^ - 20% 



a^zb 10% 


±0 


±0 


±0 


+ 


+ 


+ 








u'e ± 10% 


±0 


+ 




±0 


+ 




±0 


+ 






10.57 


9.51 


11.63 


12.68 


11.63 


13.74 


8.46 


7.40 


9.51 


M [IQ^^Mq] 


1.80 


1.98 


1.62 


1.98 


2.18 


1.79 


1.62 


1.79 


1.46 


At [years] 


9.53 


8.97 


9.98 


12.03 


11.53 


12.41 


7.21 


6.60 


7.72 



Table 4: A = 



aE± 10% 


±0 


±0 


±0 


+ 


+ 


+ 








± 10% 


±0 


+ 




±0 


+ 




±0 


+ 






11.86 


10.67 


13.05 


14.23 


13.05 


15.42 


9.49 


8.30 


10.67 


M [10^3^0] 


1.68 


1.84 


1.51 


1.84 


2.03 


1.66 


1.51 


1.66 


1.36 


At [years] 


8.70 


8.20 


9.10 


10.97 


10.53 


11.30 


6.59 


6.04 


7.05 
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4 Conclusion 



In this paper, we have computed the time delay caused by a spherical mass, a cluster of galaxies, 
in presence of a cosmological constant A using the Einstein-Straus solution, which is the appro- 
priate framework for taking into account the precession of the observer and the effect of the other 
masses of the universe in the form of a homogeneous isotropic dust, the observer being taken 
comoving with the dust. We have applied our results to the lensed quasar SDSS J1004-I-4112. 
We have computed the time delay between the images C and D of the quasar, which are the 
most aligned images with the lens, and obtained results compatible with the lower bound given 
by Fohlmeister [21j. Our predictions of the time delay range from 6 to 13 years. 

It is worthwhile to compare our results with those of previous computations performed by 
Schiicker and Zaimen [T7] and by Kawano and Oguri [22j. Schiicker and Zaimen computed the 
time delay in the framework of the Kottler solution and obtained predictions ranging from 13 
to 28 years, using however a different mass of the lens. On the other hand, Kawano and Oguri 
obtained a time delay of 10 years. 

In addition to the hypothesis of sphericity we have made in our calculations an additional 
assumption: we have supposed implicitly that the photons don't penetrate the mass distribution 
region, since we have used only the exterior Kottler solution inside the vacuole. It is then 
worthwhile to repeat the calculations taking into account that the photons can penetrate the 
interior of the mass distribution, where an interior Kottler solution must be used |23] . 
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